jji 

1 
1 

. 


re 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

NUTRIENT  STUDY  OF  MESOSCALE  THERMAL 
FEATURES  OFF  POINT  SUR,  CALIFORNIA 

by 

Walter  Elof  Hanson,  Jr. 

September  1980 

The 

sis  Advisor:           Eugene  D.  Traganza 

Approved  for  public  release;  distribution  unlimited 


T197841 


IINri.ASSTFTEn 


?"«».,, 


SECURITY   CLASSIFICATION   OF    THIS  »«fil   (W*,mt,  Dmtm  Cni»r«4) 


REPORT  DOCUMENTATION  PAGE 


AT     NUMll» 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

t.      MICl'lCNT'S  CATALOG   NUHICM 


RERO 


2.  OOVT   ACCCMIOM  NO 


4       Ti  Tl£    .  and  luBiiMti 


Nutrient  Study  of  Mesoscale  Thermal 
Features  off  Point  Sur ,  California 


••    Type  of  rsrort  *  period  covered 

Master's  Thesis; 
September  1980 


•  •  PERFORMING  org.  rerort  NUMBER 


7.   »uThO»Hi 


A.  CONTRACT  OR  SRANT  NUMICMM) 


Walter  Elof  Hanson,  Jr 


io.   rrooran  element,  rroject    task 

AREA  *  WORK  UNIT  NUDUM 


t      PERFORMING  ORGANIZATION   NAME   amO  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


It        CONTROLLING  OFFICE   NAME    UNO   AOOREtS 

Naval  Postgraduate  School 
Monterey,  California  93940 


12.  RERORT  OATE 

September  1980 


IS.     NUMBER  OF  RASES 


182 


U.     MONITORING   AGENCY   name   A    AOORESSTff  dtllmrmtl  Irmm  Contrmllln*  Oltlc*) 


IS.     SECURITY  CLASS,  (at  tnlm  ripon) 

Unclassified 


ISA.     DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


l«.     DISTRIBUTION   STATEMENT  lot  xht%  **pmr, , 

Approved  for  public  release;  distribution  unlimited. 


17.     OISTRISUTION  STATEMENT     at  th»  •fcarract  mtt—1  In  Mloem  20.  II  dttiaemtt  trmm  JtaRawfj 


IS.     SURRLEMENTARY   NOTES 


I*       KEY   WORDS  I  Continue  on  rmrtao  •'»•  II  n»c»t*mrr  jnW  id+nltfy  *r  mloek  numm+t) 


nutrients 
nitrate 
phosphate 
thermal  fronts 


chemical  fronts 

sea  surface  temperature 

satellite  infrared  imagery 

upwelling 


20      ABSTRACT  (Camtlntf  an  <••>»•«■••  •<««•  II  n«e«a««rr  *»•"  Hmmtltr  »r  »!•«*  R—R—J 

Thermal  patterns  with  the  appearance  of  cyclonic  motion 
and  sharp  thermal  fronts  frequently  are  seen  in  satellite 
IR  images  off  the  California  coast  near  Pt.  Sur.   These  thermal 
patterns  are  associated  with  distinctly  structured  upwelling 
systems.   In  nine  seasonal  cruises  since  December  1978,  nutrient 
fronts  are  strongly  correlated  with  thermal  fronts.   However, 
an  upwelling  which  is  well  defined  by  the  satellite  detected 


DD  ,  :°r7,  1473 
(Page  l) 


EDITION  OF   t  MOV  •■  IS  OBSOLETE 

S/N    0103-014- ««0I 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  TNIS  RAOE  (9hmn  Dmta  *n*fd) 


UNCLASSIFIED 


(teuwty  CL«Miy*C»T)»i|  Q»  Twit  »♦«»<•'««»—  r»..«  x.t**.* 


thermal  front,  may  not  be  nutrient  rich.   A  strong  inverse 
linear  correlation  between  nutrients  and  temperature  though 
does  exist  for  upwelling  within  an  early  stage  of  develop- 
ment.  It  is  feasible,  by  utilizing  limited  in.  situ  data, 
to  infer  nutrient  distribution  within  satellite  detected, 
surface  thermal  patterns  associated  with  upwelling. 


DD     Form       1473 

.  1  Jan  73 
S/N    0102-014-6601  tceu«<Tv  cuamiucatio*  o*  **'•  no«r»fc««  o««  ««»•»•** 


Wan  73  2  UNCLASSIFIED 

N    0102-r 


Approved  for  public  release;  distribution  unlimited 


Nutrient  Study  of  Mesoscale  Thermal 
Features  Off  Point  Sur,  California 


by 


Walter  Elof  fianson,  Jr. 

Lieutenant,  United  States  Coast  Guard 

B.S.,  United  States  Coast  Guard  Academy,  1974 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OCEANOGRAPHY 
from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  1980 


ABSTRACT 

Thermal  patterns  with  the  appearance  of  cyclonic  motion 
and  sharp  thermal  fronts  frequently  are  seen  in  satellite 
IR  images  off  the  California  coast  near  Pt.  Sur.   These 
thermal  patterns  are  associated  with  distinctly  structured 
upwelling  systems.   In  nine  seasonal  cruises  since  December 
1978,  nutrient  fronts  are  strongly  correlated  with  thermal 
fronts.   However,  an  upwelling  which  is  well  defined  by  the 
satellite  detected  thermal  front,  may  not  be  nutrient  rich. 
A  strong  inverse  linear  correlation  between  nutrients  and 
temperature  though  does  exist  for  upwelling  within  an  early 
stage  of  development.   It  is  feasible,  by  utilizing  limited 
in  situ  data,  to  infer  nutrient  distribution  within  satellite 
detected,  surface  thermal  patterns  associated  with  upwelling. 
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I.   INTRODUCTION 

This  thesis  is  a  study  of  the  linear  correlation  of  sat- 
ellite detected  sea  surface  thermal  patterns  to  nitrate  and 
phosphate  levels  in  a  coastal  upwelling  region.   In  part,  it 
considers  whether  the  chemical  mesoscale  in  upwelling  systems 
may  be  inferred  from  satellite  IR  detection  of  surface  thermal 
patterns. 

The  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on 
the  TIROS  (Television  Infrared  Observation  Satellites)  sat- 
ellite series  can  be  used  to  locate  sea  surface  temperature 
patterns  which  identify  upwelling  regions  for  in  situ  auto- 
mated biochemical  and  thermal  sampling  and  analyses.  Surface 
temperature  and  nutrient  maps  can  be  generated  using  in  situ 
measurements  and  inference  re  temperature  based  on  synoptic 
satellite  IR  imagery  and  nutrient  to  temperature  correlations. 

Previous  studies  which  used  the  same  methods  to  detect  and 
describe  the  thermal  features  (e.g.,  Traganza,  Nestor,  and 
McDonald,  1980)  showed  that  some  of  these  features  resulted 
from  upwelling  pulses  and  produced  nutrient  structures  of  the 
same  scale.   These  highly  structured  upwelling  systems  had 
sharp  thermal  and  chemical  gradients  or  fronts.   The  recurrent 
formation  of  a  "cyclonic  upwelling  system"  off  Pt.  Sur,  Cali- 
fornia [Traganza,  Conrad,  and  Breaker,  1980],  offers  a  unique 
opportunity  to  investigate  the  significance  of  frontal  systems  in 
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regions  of  upwelling.   Significantly  higher  nutrient  concen- 
trations observed  in  upwelling  systems  can  enhance  biological 
activity  in  the  proximity  of  the  fronts. 

The  ability  to  infer  the  nutrient  concentrations  associ- 
ated with  an  upwelling  system,  identified  in  satellite  IR 
imagery,  might  help  in  explaining  the  relationships  between 
the  chemical  mesoscale  and  pelagic  ecosystem  in  the  region. 
The  occurrence  of  microplanktonic  blooms  adjacent  to  sharp 
thermal  and  chemical  gradients  has  been  documented  in  many 
coastal  upwelling  regions  around  the  world  [Barber,  1980; 
Blasco,  1980;  Boyd  and  Smith,  1980;  Traganza,  Conrad,  and 
Breaker,  1980].   Because  phytoplankton  are  the  primary  pro- 
ducers in  the  pelagic  ecosystem,  their  productivity  supports 
the  foodweb  of  which  man  harvests  primarily  the  higher  level 
heterotrophs.   The  magnitudes  of  nutrient  concentrations  and 
nutrient  fluxes  in  upwelling  systems  could  provide  clues  to 
the  biological  productivity  in  such  areas.   The  associated 
thermo-chemical  gradients  or  fronts  could  become  the  foci 
for  fishing  efforts.   With  the  potential  for  biological  pro- 
duction known  and  the  growth  rate  predicted  for  an  ecosystem, 
fishing  fleets  could  be  effectively  dispatched  to  optimize 
yields . 

While  thermal  patterns  may  correlate  well  with  satellite 
observations,  correlations  for  chemical  patterns  with  IR 
imagery  are  by  inference  and  may  be  less  simple,  tending  to 
vary  as  a  function  of  the  initial  conditions,  boundary 
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conditions,  biological  processes,  and  the  dynamic  processes 
acting  over  time.   An  understanding  of  all  these  processes 
may  be  necessary  to  translate  satellite  thermal  patterns 
into  nutrient  levels. 

The  objectives  of  this  thesis  are: 

1.  To  investigate  the  linear  relationship  between  in 
situ  values  of  nutrients  (nitrate  and  phosphate)  and  temper- 
ature associated  with  coastal  upwelling  occurring  under 
varying  nutrient  and  temperature  characteristics  of  the 
ambient  water  mass  and  stages  of  development. 

2.  To  map  observed  values  of  temperature,  nitrate, 
phosphate,  and  nutrient  ratio  (nitrate/phosphate). 

3.  To  compare  the  magnitudes  and  gradients  of  tempera- 
ture and  nutrients  associated  with  the  oceanic  front  in 
different  types  of  upwelling  in  different  stages  of  develop- 
ment. 

4.  To  examine  the  relationship  between  the  surface 
temperatures  and  nutrient  concentrations  in  recently  upwelled 
water  and  to  estimate  the  depth  from  which  it  advected. 

5.  To  test  the  accuracy  of  a  satellite-aided  prediction 
in  which  the  source  water  temperature  and  nutrient  character- 
istics are  used  in  a  regression  equation  to  hindcast  surface 
nutrient  concentrations. 
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II.   METHODS 

A.  STUDY  AREA 

The  study  area  (Fig.  1)  extends  ca.  30  km  north  of  Pt. 
Sur  and  ca.  60  km  to  the  south.   It  is  bounded  by  the  122° 
30'  W  meridian  and  the  Central  California  coastline  stretch- 
ing from  Monterey  to  Cape  San  Martin. 

B.  CRUISES 

Three  cruises  were  conducted  aboard  the  Naval  Postgraduate 
School's  research  vessel  R/V  ACANIA,  27  and  28  September  1979; 
29  and  30  November  1970;  and  10  and  11  June  1980.   Positions 
were  determined  by  LORAN-C  which  has  good  coverage  in  this 
region.   Running  plots  of  temperature  and  nitrate  and  phosphate 
concentrations  were  maintained  and  used  to  ensure  coverage  of 
features  of  interest.   Particular  attention  was  given  to  the 
location  of  thermal  gradients  observed  in  satellite  images 
received  prior  to  each  cruise. 

C.  SURFACE  SEAWATER  SAMPLING  SYSTEM 

Seawater  was  continuously  sampled  from  a  depth  of  2.5m. 
The  September  and  November  cruises  pumped  seawater  from  a 
keel  intake  via  the  ship's  pump  to  the  shipboard  laboratory. 
The  June  cruise  used  an  overboard  nylon  hose  intake  at  2.5m, 
which  was  attached  to  a  deck-mounted  Wilden,  M-4,  air  oper- 
ated diaphragm  pump  (Wilden  Pump  and  Engineering  Company, 
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Fig.  1.      Bathymetry   of   the   study  area  with  grid  superimposed. 
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Colton,  California  92324)  .   An  air  debubbler  was  used  in 
line  between  the  pump  and  measuring  equipment  [Technicon 
Autoanalyzer ,  AA-II  and  Turner  III  Fluorometer;  see  Bronsink, 
1980]  . 

D.  VERTICAL  PROFILE  SAMPLING 

Seven  Nansen  casts  were  conducted  also  on  the  November 
cruise.   Seawater  samples  were  collected  at  25,  50,  100,  200, 
300,  and  400  meters  for  analysis  of  temperature,  salinity, 
and  nitrate. 

E.  TEMPERATURE 

A  thermistor  sensed  the  temperature  of  the  seawater 
entering  the  keel  intake.   A  strip  chart  recorder  continu- 
ously plotted  it.   The  equipment  was  calibrated  and  monitored 
every  30  minutes  by  bucket  thermometer  readings  of  the  sea 
surface.   The  mean  difference,  T(bucket)  -  T(thermistor) ,  of 
0.1°C  (std.  dev.  ±0.2°C)  was  determined  for  the  temperature. 
Sippican  expendable  bathythermograph  (XBT)  probes  were  re- 
leased every  20  minutes,  or  ca.  every  5  km,  during  the  Sep- 
tember cruise.   The  September  1979  and  June  1980  cruises 
dropped  XBT  probes  only  in  areas  of  major  interest.   A  spac- 
ing of  ca.  5  km  was  maintained  along  the  major  axis  of  the 
feature  and  during  front  transits.   For  calibration,  sea 
surface  temperature  was  measured  by  bucket  thermometer 
simultaneously  with  the  launch  of  each  XBT  probe.   The  mean 
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difference,  T(bucket)  -T(probe),  between  surface  measure- 
ments was  0.0°C  Cstd.  dev.  ±0.2°C). 

F.  NUTRIENTS 

Surface  (2.5m)  nutrient  concentrations  were  analyzed 
every  two  minutes  according  to  the  Technicon  Industrial  Method 
175-72-WM,  177-72-WM  [Anonymous,  1973]  and  100-70-WM  [Anon- 
ymous, 1978].   Nitrate  here  may  include  traces  of  in  situ 
nitrite,  since  the  nitrate  is  reduced  to  nitrite  before  meas- 
urement.  However,  according  to  Paulson  [1972]  there  is  little 
or  no  nitrite  at  the  surface  sampling  depth  in  this  area. 

G.  LINEAR  REGRESSION  AND  CORRELATION  COMPUTATIONS 
Regression  lines  were  generated  for  nitrate  versus 

phosphate,  nitrate  versus  temperature,  phosphate  versus  tem- 
perature, and  nutrient  ratio  (nitrate/phosphate)  versus  tem- 
perature.  The  linear  regression  analysis  used  the  NORLSQ 
library  subroutine  from  the  W.R.  Church  Computer  Center  at 
the  Naval  Postgraduate  School. 

The  population  correlation  coefficient  was  derived  by 
the  equation: 

nZx.y.  -  (Zx.)(Zy.) 
r  =  ±—± ± ± 

C[nSx.2  -  (Ix.^nnly/5  -  (ly^2])* 

When  a  "none  detected"  condition  occurred  due  to  concentra- 
tions below  the  sensitivity  of  the  instrument,  the  nutrient 
and  temperature  data  were  not  used  in  the  calculations. 
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H.   SATELLITE  IMAGERY 

The  AVHRR  has  a  temperature  resolution  of  .5°C  within 
the  IR  (10.5  ym  to  12.5  ym)  spectral  window  and  a  spatial 
resolution  of  1.1  km.   The  thermal  energy  measured  by  the 
AVHRR  is  primarily  a  function  of  radiation  emitted  from  the 
upper  millimeters  of  the  sea  surface,  clouds,  and  atmospheric 
gases.   The  greatest  errors  in  estimating  sea  surface  temper- 
ature are  induced  by  atmospheric  attenuation  of  the  thermal 
energy  due  to  absorption  and  scattering.   The  effect  of 
absorption  by  moisture  in  the  atmosphere  presents  the  most 
severe  problem  because  of  the  magnitude  of  absorption  and  the 
variability  of  the  moisture  field.   At  mid-latitudes  over  the 
ocean,  this  attenuation  may  result  in  measurements  that  are 
2  to  6°C  lower  than  actual  sea  surface  skin  temperatures 
[Maul  and  Sidran,  1973]. 

The  AVHRR  has  a  temperature  response  from  ca.  -90  to  60°C 
The  measured  thermal  energy  is  normally  displayed  as  gray 
tones  on  photographic  film  to  produce  images  interpretable 
as  clouds,  land,  and  water.   The  colder  temperatures  within 
these  images  are  assigned  lighter  shades  of  gray.   Image 
enhancement  reduced  the  range  of  temperatures  to  those  rep- 
resentative of  the  sea  surface,  while  retaining  all  the  gray 
tones  [see  Methods:  Satellite  Assisted  Sampling  Strategy]. 
In  the  enhancement  of  the  TIROS-N  images  (Plates  1,  2,  and 
3)  a  temperature  range  representative  of  the  in  situ  range 
of  values  was  assigned  to  shades  of  gray  with  white  equal 
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to  ca.  7°C  and  black  ca.  22°C.   Studying  satellite  images 
alone  is  subjective  because  of  the  inability  to  assign  tem- 
perature values  representative  of  the  sea  surface  to  each 
gray  shade.   However,  the  relative  temperature  values  and 
thermal  gradients  are  more  representative  of  in  situ  measure- 
ments within  the  same  image  than  are  absolute  values.  Because 
of  the  atmospheric  effects  on  the  thermal  energy  reaching  the 
satellite,  these  measurements  can  best  be  described  as 
"apparent  temperatures." 

I.   SATELLITE  ASSISTED  SAMPLING  STRATEGY 

Approximately  one  week  prior  to  each  scheduled  departure, 
Mr.  Breaker,  the  staff  oceanographer  at  the  National  Environ- 
mental Satellite  Service  (NESS) ,  in  Redwood  City,  California, 
would  be  alerted.   He  closely  monitored  images  of  TIROS-N 
satellite  series  for  oceanic  features  of  interest  in  the 
study  area.   The  images  were  computer  enhanced  by  him  to 
better  delineate  sea  surface  thermal  patterns.   All  images 
had  the  non-linear  distortion  associated  with  the  earth's 
curvature  and  rotation  removed.   The  day  prior  to  the  cruise 
and  sometimes  during  the  cruise,  Mr.  Breaker  sent  updated 
satellite  information  regarding  the  approximate  location, 
center,  size,  and  orientation  of  the  feature.   When  possible, 
the  most  useful  and  most  recent  images  were  delivered  prior 
to  sailing.   All  proved  invaluable  in  locating  features  and 
planning  sampling  strategy. 
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J.   SATELLITE  ASSISTED  MAPPING  OF  SURFACE  TEMPERATURE 

Two  objective  temperature  fields  were  drawn  by  the  CONISD 
computer  program;  one  field  was  formed  from  only  in  situ 
temperature,  the  other  included  data  from  satellite  imagery 
as  well. 

The  size  of  the  study  area  was  small  enough  that  the 
spherical  distortion  attributed  to  the  earth's  curvature  was 
assumed  to  be  negligible.   Therefore,  the  author  constructed 
a  grid  on  the  assumed  planar  surface  (see  Appendix  C)  with 
the  head  of  the  Sur  submarine  canyon  as  its  origin  (Fig.  1). 
Each  in  situ  data  point  had  cartesian  coordinates  assigned. 
Using  this  annotated  data,  the  CONISD  library  subroutine  from 
the  W.R.  Church  Computer  Center  drew  surface  contours.   Be- 
cause the  contouring  computations  involved  a  triangulation 
process  between  data  points  nearest  the  limits  of  the  grid 
and  the  study  area  boundaries,  it  was  necessary  to  assign 
ocean  and  near  coast  boundary  values  that  were  realistic  to 
generate  a  surface  map  which  well  represented  the  external 
limits  of  the  surface  feature.   These  values  were  inferred 
from  the  thermal  patterns  visible  in  satellite  IR  images. 
Enhanced  satellite  images  taken  during  or  nearest  the 
period  of  in  situ  observation  and  having  a  cloud-free  view 
of  the  study  area  were  used.   These  images  were  taken  to  NESS 
to  geographically  locate  the  surface  thermal  patterns  associ- 
ated with  the  upwelling  system  on  a  transverse  mercator 
projection  of  the  study  area.   The  author  outlined  shifts  in 
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the  gray  tones,  visible  to  the  naked  eye.   These  lines  rep- 
resented a  subjectively  determined  change  in  the  gray  scale. 
A  uniform  distribution  of  atmospheric  moisture  over  the  small 
area  of  the  surface  feature  was  assumed.   The  change  in  gray 
scale  could  therefore  be  related  to  synoptic  variations  in 
apparent  surface  temperature.   Mr.  Breaker,  with  his  exper- 
tise in  using  the  Baush  and  Lomb  Zoom  Transfer  Scope,  trans- 
cribed the  apparent  thermal  patterns  onto  the  transverse 
mercator  projection.   These  patterns  were  then  transferred 
by  the  author  onto  the  grid  developed  for  the  study  area. 

Using  a  surface  temperature  map  generated  only  from  in 
situ  measurements,  a  plot  of  the  ship's  track  (to  know  where 
and  when  _in  situ  measurements  were  taken)  and  the  synoptic 
thermal  patterns  inferred  from  the  satellite  image  of  the 
feature,  the  thermal  patterns  within  the  area  of  concentrated 
in  situ  measurements  were  assigned  temperature  values  equiva- 
lent to  those  on  the  surface  temperature  map.   The  satellite 
detected  and  in  situ  thermal  patterns  matched  well,  within 
±. 5°C.   Where  these  thermal  patterns  approached  the  coast, 
those  values  inferred  above  were  applied  to  the  near  coast 
boundary.   For  objective  (machine)  contouring  the  oceanic 
temperature  field  beyond  6  km  from  the  seaward  outline  of  the 
surface  thermal  feature  was  arbitrarily  assigned  a  uniform 
value.   The  oceanic  boundary  formed  a  curve  at  a  constant 
distance  (ca.  6  km)  seaward  of  the  surface  thermal  feature. 
The  warmest  in  situ  measurement  made  in  this  oceanic  region 
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was  assigned  to  the  oceanic  boundary.   With  the  inferred 
near  coast  and  ocean  boundary  values  and  in  situ  data, 
another  surface  temperature  map  was  generated  by  the  CONISD 
subroutine.   Hand  smoothing  was  applied  to  reduce  the  dis- 
tortion thought  to  be  caused  by  advection  or  other  dynamic 
processes  over  the  time  period  of  the  study.   These  temper- 
ature surface  maps  appear  in  Figures  10,  20,  and  34. 

K.   SURFACE  NUTRIENT  MAPS 

The  sea  surface  nutrient  maps  (Figs.  11,  12,  13,  21,  34, 
35,  and  36)  were  generated  in  similar  fashion  to  the  sea  sur- 
face temperature  map.   Each  in  situ  data  point  was  assigned 
cartesian  coordinates  on  the  study  area  grid.   The  near  coast 
and  ocean  boundary  values  for  temperature  were  put  into  the 
regression  equations  for  nitrate/temperature,  phosphate/ 
temperature,  and  nutrient  ratio/temperature.   The  respective 
nitrate,  phosphate,  and  nutrient  ratio  output  were  assigned 
to  the  coordinates  of  the  input  temperature.   Using  this 
annotated  In  situ  and  temperature  inferred  data,  the  CONISD 
library  subroutine  drew  surface  contours.   Hand  smoothing 
was  applied  to  reduce  the  distortion  thought  to  be  caused  by 
advection  or  other  dynamic  processes  over  the  time  period  of 
the  study. 

L.   DEFINITION  OF  FRONT  AND  CALCULATION  OF  ITS  GRADIENTS 
For  the  purpose  of  this  thesis,  the  oceanic  front  is 
defined  as  the  area  where  the  thermal  and  chemical  gradients 
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reach  their  maximum  change  per  kilometer  [see  Discussion: 
The  Relationship  Between  Source  Water  Characteristics  and 
the  Oceanic  Front] .   In  order  to  locate  the  oceanic  front 
and  assess  its  gradients,  the  following  procedures  were 
carried  out. 

The  incremental  change  per  kilometer  of  nutrient  concen- 
trations and  temperature  was  determined  from  the  in.  situ 
values  observed  along  the  ship's  track.   The  elapsed  distance 
midway  between  the  magnitudes  used  to  compute  this  spatial 
gradient  was  assigned  as  the  position  on  the  ship's  track  at 
which  the  change  occurred  (Figs.  14,  22,  and  38).   This  tech- 
nique was  used  to  take  advantage  of  the  high  spatial  resolu- 
tion (better  than  0.6  km)  along  the  ship's  track  which  showed 
microstructure  fronts  within  the  feature  as  well  as  the 
oceanic  front. 

The  magnitude  of  the  gradient  of  the  front  was  computed 
after  correcting  for  the  angle  at  which  the  ship  transited 
the  front.   It  was  assumed  that  the  front  parallels  isolines 
of  temperature  and  nutriend  concentrations  on  the  machine- 
contoured  charts.   The  angle  of  incidence  the  ship's  track 
intersected  these  isolines  at  the  vicinity  of  the  front  was 
measured  using  surface  nutrient  and  temperature  maps.   Apply- 
ing the  law  of  sines,  the  distance  along  the  normal  was  found 
and  used  to  compute  the  magnitude  of  the  gradient. 
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III.   RESULTS 

A.   SEPTEMBER  1979  CRUISE 

This  cruise  investigated  and  verified  the  existence  of  a 
satellite  sensed  wedge- like  surface  thermal  feature  extending 
100  km  SW  from  Pt.  Sur  (Plate  1).   From  satellite  images  re- 
ceived from  18  to  26  September,  it  appeared  that  this  feature 
formed  off  Pt .  Sur  on  18  September.   By  20  September  the  apex 
of  this  wedge- like  feature  had  acquired  a  northward,  anti- 
cyclonic  curvature.   Between  20  and  26  September  this  feature 
became  diffuse  ,    a  new  coastal  upwelling  event  commenced 
and  a  new  feature  was  superimposed.   A  satellite  imagery 
summary  is  presented  in  Figure  2.   The  best  image  of  the 
feature  (Plate  1)  was  taken  at  2255  GMT,  26  September,  _ca-  2 
hours  prior  to  the  ship's  departure  for  the  study  area. 
Winds,  as  measured  by  the  ship's  anemometer,  blew  from  the 
northwest  averaging  between  5  and  7  m*  sec 

The  objectives  of  this  cruise  were  to  obtain  surface 
samples  from  the  feature  and  its  oceanic  surroundings,  to 
sample  both  day  and  night  to  minimize  diel  effects,  and  to 
maximize  time  near  the  cold  core  of  the  feature.   A  seven- 
pointed  star  track  was  planned  (Fig.  3)  using  thermal  pattern 
information  derived  from  satellite  IR  imagery.   Temperature, 
nitrate,  and  phosphate  are  plotted  over  elapsed  distance  in 
Figure  4.   Correlation  coefficients  of  r  =  0.78  for  nitrate 
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Plate  1.   TIROS-N  satellite  IR  image  of  the  California  coast 
for  26  September  19  80. 
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Fig.  2.   Satellite  features  observed  from  20  September  to 
2  October  1979  [Johnson,  1980]. 
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SHIP'S    TRACK 
27-28    SEPTEMBER   1979 

SCALE        t   :    750,000 

OISTANCES    IN   KILOMETERS 

TIME    :    GMT 
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Fig.  3.   Track  of  the  September  1979  cruise  and  outline 
(dashed  line)  of  the  oceanic  front. 
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Fig.  4.   Nitrate,  phosphate,  and  sea  surface  temper 
ature  versus  elapsed  distance  along  the 
track  of  the  September  1979  cruise. 
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Fig.  4  (cont'd)   Nitrate,  phosphate,  and  sea  surface 
temperature  versus  elapsed  distance  along 
the  track  of  the  September  1979  cruise. 
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to  phosphate,  r  =  -0.91  for  nitrate  to  temperature,  and 
r  =  -0.66  for  phosphate  to  temperature  were  obtained  (Table  I). 
Nutrient  ratio  and  temperature  are  plotted  over  elapsed  dis- 
tance in  Figure  5.   The  linear  correlation  coefficient  between 
nutrient  ratio  and  temperature  is  r  =  -0.50. 

The  linear  regression  analysis  of  nitrate/phosphate, 
nitrate/temperature,  phosphate/temperature,  and  nutrient 
ratio/temperature  (Figs.  6,  7,  8,  and  9)  yielded  slopes  of 
13.85,  -2.66,  -1.00,  and  -7.72  respectively,  and  x-axis  in- 
tercepts of  0.45  yM  phosphate,  16.68°C,  16.77°C,  and  16.53°C 
respectively.   The  x- intercept  indicated  that  the  environment 
is  nitrate  limited  and  determined  the  warmest  temperature  for 
which  these  regression  lines  could  be  applied. 

The  major  axis  of  the  upwelling  system  appeared  to  parallel 
and  be  located  in  the  vicinity  of  the  Sur  submarine  canyon 
(Figs.  10,  11,  12,  and  13).   The  largest  area  of  the  feature 
and  the  part  which  appears  to  curl  cyclonically  is  south  of 
the  submarine  canyon.   The  leading  edge  of  the  curl  is  located 
between  the  Sur  and  Lucia  submarine  canyons  over  the  shelf 
slope,  beyond  the  1500  meter  depth  contour. 

The  incremental  change  per  kilometer  of  temperature  and 
nutrient  concentrations  are  plotted  over  elapsed  distance  in 
Figure  14.   The  oceanic  front  appears  to  be  located  at  elapsed 
distance  13,  18,  39,  67,  166,  200,  312,  361,  and  414  km 
(Table  II) 
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Fig.  5.  Nutrient  ratio  and  sea  surface  temperature 
versus  elapsed  distance  along  the  track  of 
the  September  1979  cruise. 
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Fig.  5  (cont'd)   Nutrient  ratio  and  sea  surface 
temperature  versus  elapsed  distance  along 
the  track  of  the  September  1979  cruise. 
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Fig.  6.   Nitrate  versus  phosphate  for  the  September  1979 
cruise. 
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Fig.  7.   Nitrate  versus  temperature  for  the  September  1979 
cruise. 
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Fig.  8.   Phosphate  versus  temperature  for  the  September  1979 
cruise. 
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Fig.  9.   Nutrient  ratio  versus  temperature  for  the 
September  1979  cruise. 
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SEA     SURFACE     TEMPERATURE 
27-28    SEPTEMBER    1979 

SCALE        1   :    750.000 

DISTANCES    IN    KILOMETERS 

CONTOURS    IN    'C 
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Fig.  10.   Sea  surface  temperature  map  for  the  September  1979 
cruise  (contour  interval,  0.5°C).   Map  generated 
from  in  situ  data  aided  by  IR  imagery. 
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NITRATE 
27  -  28    SEPTEMBER    1979 

SCALE        t  :    750.000 

DISTANCES    IN   KILOMETERS 
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Fig.  11   Surface  nitrate  map  for  the  September  1979  cruise 
(contour  interval,  1  uM  nitrate).   Map  generated 
from  in  situ  data  aided  by  inferences  from  IR 
imagery. 
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PHOSPHATE 
27-28    SEPTEMBER   1979 

SCALE        1   :    730,000 
0ISTANCE5    IN   KILOMETERS 
CONTOURS    IN    ^M 
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Fig.  12    Surface  phosphate  map  for  the  September  1979  cruise 
(contour  interval,  0.2  yM  phosphate).   Map  gener- 
ated by  in  situ  data  aided  by  inferences  from  IR 
imagery. 
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Fig.  13   Surface  nutrient  ratio  map  for  the  September  1979 
cruise  (contour  ratio  interval,  1).   Map  generated 
by  in  situ  data  aided  by  inferences  from  IR  imagery. 
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Fig.  14   Incremental  change  per  kilometer  of  temperature, 
nitrate,  and  phosphate  versus  elapsed  distance 
along  the  track  of  the  September  1979  cruise. 
Frontal  transit  occurred  at  peaks. 
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Fig.  14  (cont'd)   Incremental  change  per  kilometer  of 
temperature,  nitrate,  and  phosphate  versus 
elapsed  distance  along  the  track  of  the  September 
1979  cruise.   Frontal  transit  occurred  at  peaks. 
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The  thermal  and  chemical  gradients  are  both  sharp  and 
well  defined  despite  a  narrow  spread  of  surface  nutrient  and 
temperature  values  between  the  ocean  and  the  upwelling  system, 
0.01  to  4.36  uM  nitrate,  0.03  to  0.96  uM  phosphate,  and  17.05 
to  14.80°C.   The  oceanic  front  seemed  to  stay  in  close  prox- 
imity to  the  15.9°C,  1.7  yM  nitrate,  and  0.6  uM  phosphate 
isolines.   The  oceanic  front  is  outlined  in  Figure  3. 

Figure  15  shows  the  vertical  thermal  structure  of  the 
upwelling  system.   One  hundred  forty  XBT  probes  were  used 
with  an  average  spacing  of  4  km.   The  thermal  feature  is 
evident  from  15  to  26,  43  to  64,  197  to  246,  372  to  383,  and 
408  to  428  km.   The  thermocline  under  the  upwelling  system 
was  continuous  but  upwarped  to  within  20m  of  the  surface. 

B.   NOVEMBER  19  7  9  CRUISE 

The  staff  oceanographer  at  NESS  began  monitoring  satellite 
images  of  the  study  area  on  9  November.   He  noticed  a  discern- 
ible though  weak  thermal  feature  on  11  November.   It  extended 
35  km  southwest  of  Pt.  Sur.   From  satellite  observation,  the 
thermal  feature  appeared  to  reach  its  maximum  seaward  extent 
of  at  least  60  km  on  20  November.   By  27  November  it  had 
become  diffuse  and  by  28  November  was  indistinguishable.   Due 
to  intermittent  cloud  cover,  no  satellite  information  pertain- 
ing to  a  specific  feature  was  therefore  available  when  the 
ship  departed  28  November.   A  prior  cruise  in  the  area  from 
26  to  27  November  detected  a  frontal  feature  from  in  situ 

temperature  measurements. 
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Fig.  15  Vertical  temperature  sections  along  the  track  of 
the  September  1979  cruise  [Johnson,  1980] . 
Vertical  lines  represent  XBT  drops. 
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Fig.  15  (cont'd)   Vertical  temperature  sections  along  the 
track  of  the  September  1979  cruise  [Johnson,  1980]. 
Vertical  lines  represent  XBT  drops. 
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Fig.  15  (cont'd)  Vertical  temperature  sections  along  the 
track  of  the  September  1979  cruise  [Johnson,  1980] 
Vertical  lines  represent  XBT  drops. 
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Fig.  15  (cont'd)   Vertical  temperature  sections  along  the 
track  of  the  September  1979  cruise  [Johnson,  1980] 
Vertical  lines  represent  XBT  drops. 
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As  it  happened,  another  distinctive  thermal  pattern  was 
seen  later  in  images  on  28  November.   This  feature  extended 
west  of  Pt .  Sur  and  grew  to  maximum  size  by  6  December.   It 
advected  north  and  apparently  merged  with  upwelling  off 
Monterey  Bay  around  14  December  (Fig.  16).   Despite  inter- 
mittent cloud  cover,  the  best  satellite  image  was  taken  during 
the  study  period  at  2310  GMT,  29  November  (Plate  2).   The 
average  wind  was  3  to  5  m-sec   from  the  northwest. 

The  cruise  plan  had  the  ship  transit  south  paralleling 
the  coast.   When  the  front  was  penetrated,  an  expanding  square 
search  would  be  followed  to  delineate  the  shape,  center, 
orientation,  and  size  of  the  feature.   A  ladder  search  pattern 
was  then  to  be  executed  to  provide  additional  transits  across 
the  major  axis  of  the  feature.   When  the  expanding  square 
search  located  a  small  feature  south  of  Pt .  Sur,  the  ship 
transited  north  along  a  front  and  then  began  a  ladder  search 
seaward.   The  first  leg  transected  the  major  axis  of  a  large 
feature  located  west  of  Pt.  Sur.   After  the  general  location 
of  a  sharp  thermal  gradient  was  found  about  the  feature,  an 
expanding  square  search  was  begun  in  the  vicinity  of  the 
leading  edge  (Fig.  17). 

The  data  herein, presented  is  associated  with  those  meas- 
urements taken  during  this  ladder  search  and  final  expanding 
square.   In  addition,  seven  oceanographic  stations  were 
occupied  near  the  end  of  the  study  period.   Nansen  casts  were 
made  to  sample  temperature,  salinity,  and  nitrate. 
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Fig.  16   Satellite  feature  observed  from  19  November  to 
14  December  1979  [Johnson,  1980]. 
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Plate  2.   TIROS-N  satellite  image  of  the  California 
coast  for  29  November  1979. 
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Fig.  17   Track  of  the  November  1979  cruise  and  outline 
(dashed  line)  of  the  oceanic  front. 
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Because  of  contamination  in  the  combined  working  reagent, 
no  reliable  phosphate  concentrations  could  be  measured  on 
the  Autoanalyzer . 

The  inverse  correlation  between  nitrate  and  temperature 
[Fig.  18  and  19)  was  strong,  r  =  -0.93.   The  regression 
line  of  nitrate  to  temperature  had  a  slope  of  -3.24  and  an 
x-intercept  of  14.74°C.   The  feature  extended  seaward  in  a 
southwesterly  direction  from  a  point  near  the  coast  ca-  20  km 
north  of  Pt.  Sur.   At  a  position  ca«  20  km  west  of  Pt.  Sur, 
it  appeared  to  curl  cyclonically  (Figs.  20  and  21).   Unlike 
the  September  feature,  this  upwelling  system  did  not  appear 
to  extend  beyond  the  shelf  slope  to  depths  greater  than  1200m 

The  incremental  change  per  kilometer  of  temperature  and 
nitrate  are  plotted  over  elapsed  distance  in  Figure  22.   The 
oceanic  front  appears  to  be  transited  at  elapsed  distances 
151,  196,  203,  231,  257,  and  270  km. 

The  thermal  and  nitrate  gradients  are  sharp  and  well 
defined.   The  steepest  appear  to  occur  near  the  equatorward 
edge  of  the  curl.   Compared  with  the  September  cruise,  the 
nitrate  gradient  at  the  front  is  greater  for  the  December 
feature,  whereas  the  temperature  gradients  show  no  great 
variation  (Table  II) .   The  spread  of  the  surface  temperature 
and  nitrate  values  between  the  ocean  and  the  upwelling  fea- 
ture, 15.50  to  11.85°C  and  0.12  to  9.92  yM  nitrate,  was 
greater  in  magnitude  than  the  September  cruise.   The  oceanic 
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Fig.  18   Nitrate  and  sea  surface  temperature  versus  elapsed 

distance  along  the  track  of  the  November  1979  cruise. 
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Fig.  19  Nitrate  versus  temperature  for  the  November  1979 
cruise. 
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Fig.  20   Sea  surface  temperature  map  for  the  November  1979 
cruise  [contour  interval  0.5°C).   Map  generated 
from  in  situ  data  aided  by  IR  imagery. 
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Fig.  21   Surface  nitrate  map  for  the  November  1979  cruise 
(contour  interval,  1  yM  nitrate).   Map  generated 
from  in  situ  data  aided  by  inferences  from  IR 
imagery. 
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Incremental  change  per  kilometer  of  temperature  and 
nitrate  versus  elapsed  distance  along  the  track  of 
the  November  1979  cruise.   Frontal  transit  occurred 
at  peaks. 
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front  seemed  to  closely  approximate  the  14.0°C  and  2.5  uM 
nitrate  isolines.   The  oceanic  front  is  noted  in  Figure  17. 

Figure  23  shows  the  vertical  profile  of  the  upwelled 
feature  based  on  40  XBT  profiles.   The  vertical  thermal 
structure  is  apparent  from  elapsed  distance  179  to  216,  226 
to  247,  275  to  290,  and  298  to  319  km.   The  thermocline  under 
the  upwelling  system  is  continuous  and  shows  no  great  pertur- 
bations . 

The  data  obtained  from  the  seven  Nansen  casts  appears 
in  Figures  24  and  25.   Because  of  the  increase  in  nitrite  as 
depth  increases,  the  actual  concentration  of  nitrate  below 
the  surface  layer  can  be  less  than  the  values  recorded  (see 
Methods:  Nutrients). 

Station  7  was  taken  near  the  major  axis  of  the  upwelling 
system,  whereas  station  2  sampled  the  warmer  oceanic  water 
to  the  south  of  the  feature.   Stations  3,  4,  5,  and  6  were 
all  taken  in  close  proximity  to  the  oceanic  front.   An  in- 
trusion of  warmer  oceanic  water  into  the  central  region  of 
the  feature  appeared  to  have  been  recorded  at  station  1.   The 
Nansen  cast  data  showed  that  the  upwelling  system  apparently 
corresponded  to  a  26.4  a      tongue  of  water  upwarping  close  to 
the  coast  and  travelled  seaward  near  the  surface.   The  nitrate 
concentrations  seem  to  closely  agree  with  the  density  struc- 
ture. 
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Fig.  23  Vertical  temperature  sections  along  the  track  of 

the  November  1979  cruise  [Johnson,  1980] .   Vertical 
lines  represent  XBT  drops;  heavy  lines  represent 
bottom  topography. 
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Fig.  23  (cont'd)   Vertical  temperature  sections  along  the 
track  of  the  November  1979  cruise  [Johnson,  1980]. 
Vertical  lines  represent  XBT  drops;  heavy  lines 
represent  bottom  topography. 
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Fig.  23  (cont'd)   Vertical  temperature  sections  along 

the  track  of  the  November  1979  cruise   [Johnson, 
1980].   Vertical  lines  represent  XBT  drops. 
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Fig.  25   T  -  S  and  T  -  N  diagrams  of  seven  stations  on  the 
November  1979  cruise. 
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C.   JUNE  1980  CRUISE 

Satellite  monitoring  of  the  waters  off  Pt .  Sur  on  9  June 
indicated  a  pulse  of  upwelled  water  had  surfaced  just  before 
the  ship  departed  Monterey  en  route  to  the  study  area. 
Plate  3  was  the  best  image  of  the  surface  thermal  feature, 
taken  at  2242  GMT  9  June  1980,  four  hours  prior  to  the  com- 
mencement of  the  study.   The  wind  during  the  study  period, 
June  10  and  11,  averaged  between  11  and  13  m*sec   and  was 
blowing  from  the  northwest. 

With  the  feature  located  by  satellite  thermal  patterns, 
a  five-pointed  star  track  was  planned  so  that  crossings  of 
the  feature  and  coastal  upwelling  could  be  made.   By  halving 
the  ship's  speed  during  front  transits,  finer  resolution  data 
was  to  be  obtained  to  compare  the  poleward  and  equatorward 
fronts  for  differences  in  thermal  and  nutrient  gradients. 
The  cruise  plan  was  modified  because  of  rough  seas.   After 
three  crossings  of  the  feature,  the  front  analysis  was  given 
priority.   A  modified  ladder  track  maneuvered  the  ship  back 
and  forth  across  the  oceanic  front.   Before  concluding  the 
study,  the  ship  transected  the  upwelling  along  its  major  axis. 
The  ship's  track  is  shown  in  Figure  26. 

The  correlation  coefficients  were  r  =  0.96  for  nitrate 
to  phosphate,  r  =  -0.96  for  nitrate  to  temperature,  r  =  -0.92 
for  phosphate  to  temperature,  and  r  =  -0.96  for  nutrient  ratio 
to  temperature  (Figs.  27  and  28).   Linear  regression  analyses 
between  the  above  parameters  yielded  slopes  of  18.05,  -4.24, 
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Plate  3.   TIROS-N  satellite  image  of  the  California  coast 
for  9  June  1980. 
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Fig.  26  Track  of  the  June  1980  cruise  and  outline  (dashed 
line)  of  the  oceanic  front. 
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Fig.  27   Nitrate,  phosphate,  and  sea  surface  temperature 
versus  elapsed  distance  along  the  track  of  the 
June  1980  cruise. 
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Fig.  28  Nutrient  ratio  and  sea  surface  temperature  versus 
elapsed  distance  along  the  track  of  the  June  1980 
cruise. 
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-0.64,  and  -2.07  respectively,  and  x-axis  intercepts  of 

1.07  yM  phosphate,  13.47°C,  13.88°C,  and  13.97°C  respectively 

(Figs.  29,  30,  31,  and  32) . 

The  major  axis  of  the  thermal  and  nutrient  patterns  are 
parallel  to  and  directly  over  the  axis  of  the  Sur  submarine 
canyon  (Figs.  1,  33,  34,  35,  and  36).   The  leading  edge  of 
the  thermal  feature  appeared  to  curl  cyclonically  in  the 
vicinity  of  the  head  of  Sur  submarine  canyon.   The  crossing 
of  the  oceanic  front  appeared  to  have  occurred  at  elapsed 
distances  14,  36,  112,  158,  172,  206,  and  222  km  (Fig.  3  7). 

The  sharp  and  well  defined  thermal  and  chemical  gradients 
at  the  front  seem  to  be  less  than  in  the  November  feature  but 
greater  than  in  September  (Table  II).   However,  the  spread  of 
temperature  and  nutrient  concentrations  is  the  widest  for  the 
June  cruise;  13.77  to  9.04°C,  0.29  to  19.98  yM  nitrate,  and 
0.91  to  2.31  yM  phosphate.   The  oceanic  front  appeared  to 
closely  approximate  the  11.5°C,  8.8  yM  nitrate,  and  1.6  yM 
phosphate  isolines.   Figure  26  shows  the  oceanic  front. 

The  vertical  temperature  profile  was  generated  from 
23  XBT  profiles.   The  thermocline  under  the  upwelling  was 
discontinuous  and  appeared  to  surface  in  the  vicinity  of  the 
oceanic  front. 
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Fig.  29  Nitrate  versus  phosphate  for  the  June  1980  cruise 
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Fig.*  30  Nitrate  versus  temperature  for  the  June  1980  cruise. 
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Fig.  31   Phosphate  versus  temperature  for  the  June  1980  cruise. 
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Fig.  32   Nutrient  ratio  versus  temperature  for  the  June 
1980  cruise. 
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Fig.  33   Sea  surface  temperature  map  for  the  June  1980 

cruise  (contour  interval,  0.5°C).   Map  generated 
from  in  situ  data  aided  by  IR  imagery. 
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Fig.  34   Surface  nitrate  map  for  the  June  1980  cruise 
(contour  interval,  1  yM  nitrate).   Map 
generated  by  in  situ  data  aided  by  inferences 
from  IR  imagery. 
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Fig.  35   Surface  phosphate  map  for  the  June  1980  cruise 
Ccontour  interval,  0 .  2  jjM  phosphate).   Map 
generated  from  in  situ  data  aided  by  infer- 
ences from  IR  imagery. 
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Fig.  36   Surface  nutrient  ratio  map  for  the  June  1980  cruise 
(contour  ratio  interval,  1).   Map  generated  from 
in  situ  data  aided  by  IR  imagery. 
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TABLE  II 
GRADIENTS  AND  OBSERVED  VALUES  AT  THE  OCEANIC  FRONT 


CKIJISC 


I'O',  I . 
TIOIJ 


0  RADII:  NT 
AT    PKOIIT 


iM   NO 


/** 


H    PO„ 


km 


OBSERVED   VALUES 

AT    ri'.ONT 


T   °C     M"    "0"3  fXM   P04- 


poniTion  on  n'ntJT 
where  '-.i-aiji i.trr: ;  and 

value:,  mealuked 


DEC    79 


0.17 
0.  38 


1.09 
2.b7 


0.07 
0.16 


11.? 
11.1 


9.0 
10.0 


1.0 
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u;a«/*i 


APR   79 


0.3S 
0.34 
0.31 
0.34 
0.29 


2.13 
1.74 
1.43 
1.93 
2.33 


0.15 
0.10 
0.09 
0.13 
0.13 


10 
11 
10 
10 

11 


14, 

a, 
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11, 
12, 


1.4 
1.2 
1.5 
1.3 
1.3 


u)AR/vi 


JUN   79 


0.45 
0  .44 


2.05 
2.56 


0.12 
0.19 


13.2 
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7.1 
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0.7 
0.4 


AUG   79 


0.4.1 
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1.56 
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0.08 
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0.09 
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IV.   DISCUSSION 

A.   SATELLITE  DETECTED  THERMAL  PATTERNS  COMPARED 
TO  IN  SITU  SURFACE  MAPS 

Inferring  the  distribution  of  surface  nutrient  concen- 
trations from  surface  thermal  patterns  detected  by  satellites 
appear  feasible  from  the  correlation  of  in  situ  nutrients  and 
temperature.   In  this  study  sea  surface  nutrient  maps  were 
generated  for  comparison  with  satellite  detected  thermal 
patterns  formed  by  upwelling  systems  off  Pt.  Sur,  California. 
For  example,  the  surface  nitrate  and  phosphate  maps  for  the 
June  1980  feature  (Figs.  34  and  35)  appear  similar  to  the 
thermal  distribution  (Fig.  33)  with  regard  to  the  general 
location  of  the  nutrient  feature,  its  center,  and  the  orien- 
tation of  the  major  axis.   However,  the  nutrients  exhibited 
a  highly  structured  pattern  with  sharp  gradients  within  the 
feature  which  was  not  predictable  from  the  satellite  image. 
The  magnitude  of  the  difference  between  the  distribution 
implied  by  the  image  and  the  actual  surface  distribution 
of  nutrients  within  the  feature  varied  among  the  April,  Sep- 
tember, November  1979  and  June  1980  upwelling  systems.   (The 
other  features  lacked  sufficient  data  to  construct  surface 
nutrient  maps.)   The  thermal  and  spatial  resolution  of  the 
AVHRR  (1°C  and  1.1  km  respectively)  could  account  for  some 
of  the  difference.   However,  the  relative  ranges  of  the 
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nutrient  values  compared  to  the  range  of  the  temperature 
varied  too  much  in  different  features  at  different  times  to 
be  satisfactorily  explained  by  the  lack  of  resolution  alone. 

B.   NUTRIENT- TEMPERATURE  CORRELATION  AND 
"ACTIVE  UPWELLING"  SYSTEMS 

To  better  evaluate  these  differences,  regression  analyses 
were  applied  to  in  situ  values  of  nutrients  and  temperature 
within  the  different  features  which  formed  at  different  times 
of  the  year  (Fig.  38).   Although  there  were  variations,  the 
slopes  for  all  features  were  negative.   A  strong  inverse 
linear  correlation  existed  in  five  of  the  seven  features 
observed  since  December  1978  (Table  I) .   This  implied  that 
regression  lines  between  nutrients  and  temperature  might  be 
useful  in  describing  the  distribution  of  nutrients  with  re- 
spect to  surface  temperature  patterns.   The  variations  in 
slope  of  these  regression  lines  suggests  that  the  relation- 
ship of  nutrients  to  temperature  might  depend  on  numerous 
factors.   These  may  include  the  temperature  and  nutrient 
characteristics  of  the  source  water,  the  depth  from  which 
the  water  was  advected,  wind  stress  and  bathymetric  effects, 
surface  dynamic  processes  (viz.,  mixing,  advection,  and  heat 
exchange)  and  biological  processes.   Because  of  the  complexity 
of  the  problem,  trends  were  investigated  with  a  view  toward 
developing  indices  for  prediction  of  nutrient  distribution. 

The  strong  inverse  correlations  between  temperature  and 
nutrients  occurring  in  five  cruises  suggested  that  ' active 
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Fig.  38   Summary  of  regression  lines  for  nitrate,  phosphate, 
and  nutrient  ratio  versus  temperature  7  DEC  1978  - 
10  JUN  1980.   20  JAN,  oceanic  water;  7  DEC,  30  APR, 
8  AUG,  and  10  JUN  upwelled  water  from  below  thermo- 
cline;  remaining  cruises  "upwelled"  water  from 
above  thermocline. 
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upwelling'  systems  had  occurred.   This  was  readily  apparent 
for  the  April  1979  and  June  1980  cruises  for  which  the  ver- 
tical temperature  profiles  showed  the  thermocline  surfacing 
near  the  edge  of  the  surface  thermal  feature.   This  break, 
through  the  thermocline,  would  have  permitted  the  advection 
of  nutrient-rich,  "biochemically  new"  water  [Traganza  e_t  al . , 
1980],  to  the  surface  waters  within  the  feature.   This  seems 
to  be  supported  by  the  high  nutrient  concentrations,  high 
nutrient  ratios,  and  low  temperatures  below  10°C  observed 
inside  the  feature  (Figs.  27  and  28,  respectively). 

Although  satellite  history  of  the  November  1979  cruise 
suggests  the  observed  feature  was  sampled  during  a  stage  of 
'active  upwelling',  XBT  data  indicate  that  the  system  appeared 
not  to  have  penetrated  the  thermocline.   The  thermocline  under 
the  feature  was  continuous  and  did  not  appear  to  be  perturbed 
by  the  shallow  (ca.  50m)  upwelling  system.   This  would  have 
accounted  for  the  moderate  level  of  nitrate  observed  at  the 
surface  inside  the  feature. 

The  December  1978  and  June  1979  cruises  did  not  collect 
vertical  temperature  data.   From  the  wind  data  and  the  mod- 
erate levels  of  nitrate  and  phosphate,  the  upwelling  system 
observed  on  the  June  1979  cruise  was  assumed  ' active  ' 
[Conrad,  1980].   Because  the  nutrient  levels  were  similar 
in  magnitude  to  the  November  1979  cruise,  this  feature  could 
be  representative  of  a  shallow  upwelling. 
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As  the  December  1978  cruise  had  low  temperatures,  high 
nutrient  concentrations  and  nutrient  ratios  similar  to  those 
of  the  April  1979  and  June  1980  cruises,  the  feature  observed 
could  have  been  associated  with  a  discontinuity  in  the  ther- 
mocline  permitting  deep  upwelling.   Thus,  strong  linear 
correlations  between  nutrients  and  temperature  may  be  assumed 
for  'active  upwelling'  systems. 

C.   THE  RELATIONSHIP  BETWEEN  SOURCE  WATER 
CHARACTERISTICS  AND  THE  OCEANIC  FRONT 

The  ability  to  determine  the  initial  condition  of  the 
upwelled  water  was  considered  a  major  factor  in  prediction. 
The  nutrients  in  this  subsurface  source  water  would  provide 
the  maximum  nutrient  concentrations  observed  in  the  upwelled 
surface  water.   The  nutrient  concentrations  and  temperature 
were  expected  to  be  a  function  of  depth  and  ambient  water 
mass.   The  depth  of  the  source  water  is  most  likely  tied  to 
the  duration  and  strength  of  wind  or  other  forcing  functions, 
e.g.,  internal  waves  [Johnson,  1980]. 

Within  an  ' active  upwelling,  T  the  surface  is  replenished 
by  vertical  advection  of  source  water.   The  properties  char- 
acteristic of  the  source  water  are  confined  to  the  surface 
feature  by  the  oceanic  front.   For  every  feature  studied,  the 
maximum  gradient  analysis  of  in  situ  temperature  and  nutrients 
(see  Methods:  Definition  of  Front  and  Calculation  of  its  Grad- 
ients) showed  the  nutrient  fronts  coincided  within  +.6  km  to 
the  thermal  front  (Figs.  14,  22,  and  37).   Thus,  the  thermal 
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front  detected  in  satellite  imagery  closely  approximates 
the  nutrient  front. 

Since  nutrients  are  linearly  related  to  temperature, 
thermal  gradients  are  well  represented  within  a  satellite 
image,  and  the  nutrient  fronts  and  thermal  front  coincide, 
the  features  were  examined  to  see  if  the  sharpness  of  the 
thermal  front  might  reflect  the  range  of  nutrients  within 
the  feature. 

The  sharpness  of  the  thermal  front  was  not  found  to  be 
a  dependable  measure  of  the  nutrient  distribution  within  a 
feature.   The  thermal  gradients  at  the  fronts  of  the  December 
1978,  April  and  September  1979,  and  June  1980  features  had 
similar  magnitude,  0.12°C/km  to  0.38°C/km  (Table  II).   Within 
the  December,  April,  and  June  features  high  concentrations, 
8.3  to  23.1  yM  nitrate  and  1.0  to  2.3  yM  phosphate  were  ob- 
served.  These  surface  values  were  considered  the  result  of 
'active'  deep  upwelling.   However,  the  September  feature  had 
significantly  different  nutrient  concentrations,  0.6  to  4.3 
yM  nitrate  and  0.5  to  0.9  yM  phosphate.   The  vertical  temper- 
ature profile  of  the  September  cruise  showed  a  continuous 
though  greatly  perturbed  thermocline  under  the  surface  feature 
(Fig.  15).   The  low  nutrient  concentrations  and  nutrient 
ratios,  together  with  the  vertical  temperature  date  [see 
Johnson,  1980] ,  suggest  that  a  shallow  upwelling  or  mixing 
in  the  upper  layer  caused  the  feature.   The  satellite  history 
of  the  feature  revealed  the  thermal  patterns  were  becoming 
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indistinguishable  at  the  time  of  the  cruise.   The  lack  of 
strong  correlation  in  the  nitrate  to  phosphate  and  nutrient 
to  temperature  regressions  also  suggested  that  the  feature 
was  dissipating.   The  June,  August,  and  November  1979  fea- 
tures also  showed  differences  between  nutrient  distributions 
despite  their  similar  thermal  fronts. 

D.  NUTRIENT  DISTRIBUTION  PREDICTION   REQUIRES  GROUND  TRUTH 
Because  of  the  inability  to  describe  adequately  and  there- 
fore predict  the  nutrient  distribution  from  satellite  detected 
thermal  patterns  alone,  ground  truth  data  are  needed.   In 
that  the  linear  correlation  was  strong  between  nutrients  and 
temperature,  an  indicator  of  active  upwelling,  for  the  five 
features  mentioned  earlier,  the  sum  effects  of  the  initial 
conditions,  boundary  conditions,  and  dynamic  processes  appear 
confined  to  altering  the  slope  of  the  nutrient  to  temperature 
regression  line  and  its  intercept. 

E.  THE  FEASIBILITY  OF  SATELLITE  PREDICTION 
AIDED  BY  LIMITED  IN  SITU  SAMPLING 

A  test,  therefore,  was  conducted  to  determine  how  good  a 

satellite  prediction  could  be,  using  source  water  temperature 

and  nutrient  characteristics  in  a  regression  equation  to 

forecast  surface  nutrient  concentrations.   The  regression 

line  is  bounded  by  the  coldest  temperature  in  the  feature 

and  the  warmest  oceanic  water  temperature. 
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1.  Use  of  the  Oceanic  Water  Temperature  as  an 
Intercept  for  a  Predictive  Regression  Line 

The  temperature  (T)    intercepts  computed  from  in  situ 
nutrient/temperature  regression  equations  were  compared  with 
the  warmest  oceanic  temperature  to  determine  if  a  trend  could 
be  seen  that  would  be  useful  in  developing  a  regression  equa- 
tion to  forecast  nutrient  concentrations.   The  phosphate/ 
temperature  T- intercept  was  always  greater  (warmer)  than  the 
nitrate/temperature  T-intercept.   Because  the  oceanic  en- 
vironment appeared  to  be  nitrate-limited  (nitrate  goes  to 
zero  while  phosphate  is  still  present,  see  Figs.  6  and  29), 
the  phosphate/temperature  line  would  have  to  be  extrapolated 
to  reach  the  temperature  axis,  thus  greater  (warmer)  temper- 
atures.  This  would  be  an  over  estimate  of  the  actual  oceanic 

water  temperature.   The  difference,  T      -T.       . ,  for 
r  '   ocean   intercept' 

active  upwelling  systems  as  defined  by  nutrient  values  and 
concentrations  ranged  from  -0.21  to  +0.61°C  for  nitrate/ 
temperature  and  -0.18  to  +0.27°C  for  phosphate/temperature 
regression  equations  (Table  I) .   Within  the  thermal  resolu- 
tion (.5°C)  of  the  satellite  sensor,  the  T-intercept  might 
therefore  be  considered  equivalent  to  the  temperature  of  the 
oceanic  water. 

2.  Hypothesis  That  Nutrient/Temperature 
Relation  within  Source  Water  Is  Linear 

Strong  linear  correlations  were  seen  between  both 

nutrients  and  temperature  and  their  gradients  in  the  November 
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1979  shallow  upwelling  system  and  June  1980  deep  upwelling 
system.   From  satellite  history,  these  features  appear  to 
have  been  observed  during  an  initial  stage  of  development 
indicating  the  water  upwelled  from  its  source  depth  had  been 
in  contact  with  the  surface  for  a  day  or  less.   Therefore, 
the  source  water  nutrient  and  temperature  characteristics 
should  be  linearly  related.   If  this  is  the  case,  and  the 
source  water  regression  line  is  known  (from  recent  or  his- 
torical data),  this  relation  could  be  used,  along  with  the 
oceanic  water  temperature  as  a  T-intercept,  to  provide  a 
first  approximation  of  the  nutrient  distribution  reflected 
by  the  surface  thermal  patterns. 

In  order  to  test  this  hypothesis,  source  water  char- 
acteristics were  examined  using  T  -  S  and  T  -  N  diagrams  from 
the  shallow  upwelling  November  1979  Nansen  casts  (Fig.  25). 
The  2  and  25m  T  -  N  samples  for  station  7  (located  near  the 
coldest  water  in  the  feature)  and  stations  3,  4,  5,  and  6 
(located  near  the  front)  were  plotted  on  the  T-N  diagram  for 
station  2,  which  was  representative  of  the  oceanic  environment 
(Fig.  39).   The  25m  samples  were  tightly  clustered  on  the 
station  2  curve  between  the  25  and  50m  depths,  suggesting 
that  the  water  had  upwelled  from  that  depth.   The  2m  samples 
of  the  ocean  front  stations  were  skewed  off  the  station  2 
reference  curve  toward  warmer,  lower  nutrient  water.   However, 
station  7  varied  little  between  the  surface  and  25m  depth. 
Apparently  station  7  was  nearest  the  area  where  the  source 
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Fig.  39  T-  S  and  T-N  diagram  for  station  2;  ocean  reference, 
with  25m  (solid  circles)  and  2m  (open  circles)  data 
from  stations  3,  4,  5,  6,  and  7  added  to  identify 
source  water. 
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water  first  came  in  contact  with  the  surface.   The  ocean 
front  station  data  were  obtained  in  this  source  water  after 
it  had  been  affected  by  mixing,  heat  exchange,  and  other 
processes,  dynamical  and  biological. 

3.  The  Effects  of  Wind  Mixing  and  Heat  Exchange  to 
the  Atmosphere  on  Surface  Water  near  Oceanic  Front 

The  difference  in  2  and  25m  temperatures  could  be 

associated  in  part  with  the  mixing  caused  by  the  3  to 

5m#sec   winds  during  the  November  cruise  and  the  atmospheric 

thermal  exchange  associated  with  a  low  percentage  of  cloud 

cover  and  ca.  10  hours  of  daylight  observed.   Over  the  one 

day,  the  upwelled  water  known  to  have  been  in  contact  with 

the  surface  the  longest  (which  was  assumed  to  be  the  water 

at  the  oceanic  front)  was  predicted  to  be  ca..  0.8°C  warmer 

[James,  1966].   This  agreed  with  what  was  observed  on  the 

T  -  N  diagram  (Fig.  39) . 

4.  Determination  of  Depth  From  Which 
Source  Water  Was  Advected 

It  was  postulated  that  a  nutrient  source  depth  range 
could  be  determined  by  the  depths  at  which  the  cold  tempera- 
ture in  the  center  of  the  feature,  T  ,  and  the  temperature 
of  the  ocean  front,  T r ,  were  found  on  a  temperature  versus 
de.pth  profile,  obtained,  e.g.,  with  an  XBT.   From  these  depths 
a  range  of  nutrient  values  could  be  measured  by  .in  situ 
sampling.   The  change  in  nutrients  versus  the  change  in 
temperature  over  the  depth  range,  corresponding  to  T  -  TV , 


93 


could  be  a  first  approximation  of  the  slope,  m,  for  a  regres- 
sion equation  to  predict  nutrient  concentrations  within  the 
surface  feature.   This  slope  approximation  could  be  combined 
with  the  oceanic  water  surface  temperature  as  described 
earlier  to  compute  the  nutrient  CN)  intercept,  b,  whereby 
the  total  regression  equation,  N  =  mT  +  b,  can  be  formulated. 
The  conditions  for  which  this  approximation  would  be  best 
suited  are  those  in  which  there  were  not  thermal  exchanges 
with  the  environment  and  changes  in  nutrient  concentrations 
due  to  biological  activity.   Accuracy  of  the  application  to 
IR  satellite  data  would  depend  on  the  extent  of  agreement 
between  the  surface  temperature  inferred  from  the  satellite 
image  and  the  actual  surface  temperature. 

5.   A  Hindcast  of  Nitrate'  Concentration  Using 
Data  from  the  November  1979  Cruise 

The  November  1979  cruise  station  2  vertical  tempera- 
ture profile  showed  the  nitrate  source  depth  range  for 
11.85°C,  T  ,  and  14.0°C,  Tf ,  to  be  40  to  14m  respectively. 
The  nitrate  range  between  these  depths  was  11.7  to  3.9  yM. 
The  slope  for  this  line,  therefore,  was  ca.  -3.7.   A  second 
approximation  correcting  for  wind  mixing  and  thermal  exchange 
with  the  atmosphere  [James,  1966]  used  13.2°C  for  Tf  to  gen- 
erate a  slope  of  ca.  -3.1.   The  oceanic  water  surface  temper- 
ature used  to  predict  the  N-intercept,  b  =  -mT,  was  15.50°C. 
From  these  few  measurements,  a  predictive  regression  equation 
was  formulated. 
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The  predictive  regression  equation  computed  from  333 
in  situ  measurements  was  compared  with  the  equations  developed 
using  the  slope  derived  from  source  water  characteristics  and 
the  temperature  of  the  ocean  water  used  as  the  T-intercept. 
The  nitrate  concentrations  for  11.85°  and  14.0°C  were  computed. 

The  best  fit  linear  regression  (correlation,  r=-0.93) 
based  on  underway  in.  situ  measurements  was: 

N  =  -3.24  T  +  47.74. 

where  N  =  yM  nitrate  and  T  =  °C  temperature.   The  predictive 
regression  equation,  assuming  no  thermal  exchange  with  the 
environment,  was: 

N1  =  -3.7T  +  57.35 

The  second  approximation  which  corrected  for  wind  mixing  and 
atmospheric  thermal  exchange  was: 

N2  =  -3.1 T  +  48.05 

a.   Comparison  of  Predicted  Nitrate  Concentration 
to  That  Computed  by  Regression  Analysis  of 
Observed  Data 

The  solution  to  the  equations  are  as  follows: 

for  Tc  of  11.85°C,  N  =  9.82  yM,  H      =  13.5  yM,  and  N2  =  11.3  yM 

nitrate;  for  Tf  of  14.0°C,  N  =  2.94  yM,  N,  =  5.5  yM,  and 
1  L  N,  -  N 

N?  =  4.7  yM  nitrate.   The  percentage  of  error  ( x  100 

£        N2  -  N  N 

and  x  100) ,  associated  in  predicting  the  nitrate 

N 

concentration  at  the  cold  center  or  most  recently  upwelled 
water  appears  to  be  ca.   +37%  for  the  first  approximation, 
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N, ,  and  ca.  +15%  for  N2  which  accounted  for  wind  mixing  and 
thermal  exchange.   In  the  vicinity  of  the  ocean  front,  the 
percentage  of  error  increased;  for  TV  of  14.0°C,  N,  had  an 
error  of  c_a.  +98%  whereas  N~  had  an  error  of  ca.  +581.   In 
both  hindcasts  the  nutrient  concentrations  predicted  were 
higher  than  those  computed  from  the  equation  based  on  in  situ 
surface  measurements.   Both  hindcasts  also  deteriorated  in 
reliability  proceeding  away  from  the  cold  center  of  the  fea- 
ture.  The  predictive  equation  which  corrected  for  wind  mixing 
and  thermal  exchange  with  the  atmosphere  had  significantly 
less  error  than  the  equation  that  did  not.   Because  nitrate 
and  temperature  are  nonconservative  within  the  upper  layer 
of  the  ocean,  the  15%  error  derived  from  linear  regression 
should  be  considered  reasonable.   In  all,  the  longer  the  up- 
welled  water  is  in  contact  with  the  surface  layer  of  the 
ocean,  the  more  complex  the  interactions  with  the  physical 
and  biological  environment  may  become,  and  the  more  unreal- 
istic the  simple  predictive  equations. 

b.   Summary  of  Conditions  for  Which  This 
Test  Appeared  Satisfactory 

This  hindcast  was  applied  only  to  a  shallow  up- 
welling  system  which  satellite  history  showed  to  be  in  an 
initial  stage  of  development.   The  prediction  was  based  on 
knowledge  of  the  temperature  of  the  sea  surface  within  the 
coldest  thermal  pattern  detected,  surface  temperature  of  the 
oceanic  front  detected  by  thermal  patterns,  a  vertical 
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temperature  and  nutrient  profile  within  the  upper  200m  of 
the  ocean  seaward  of  the  oceanic  front,  and  atmospheric 
parameters  to  hindcast  ocean-atmosphere  thermal  exchange 
[James,  1966].   In  this  special  case,  given  this  limited 
ground  truth,  it  was  possible  to  hindcast  within  15%  error 
the  maximum  nutrient  concentration  in  an  upwelling  system 
detected  by  satellite  imagery.   It  therefore  seems  reasonable 
to  assume  that  at  least  the  major  patterns  of  nutrient  dis- 
tributions can,  with  further  research,  be  inferred  using 
satellite  imagery  and  limited  ground  truth  (i.e.,  data  from 
buoys  and  AXBT's) . 
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V.   CONCLUSIONS 

1.  Inference  of  nutrient  distribution  by  satellite 
detected  upwelling  systems  is  feasible. 

2.  Active  upwelling  systems  are  expected  to  have  strong 
inverse  linear  correlations  between  nutrients  and  temperature. 

3.  The  nutrient  front  position  can  be  approximated 
closely  by  the  thermal  oceanic  front. 

4.  The  nutrient  distribution  within  a  feature  can  not 
be  related  to  the  sharpness  of  the  thermal  front. 

5.  To  predict  nutrient  distributions,  ground  truth  as 
well  as  satellite  detected  thermal  patterns  are  required. 

6.  A  linear  regression  can  be  used  to  forecast  nutrient 
maxima  for  upwelling  systems  in  the  initial  stage  of  develop- 
ment aided  by  only  limited  in.  situ  data. 

7.  The  approximation  of  nutrient  concentrations  by  linear 
regression  can  be  improved  by  estimating  the  effects  of  wind 
mixing  and  thermal  exchange  with  the  atmosphere. 

8.  With  greater  knowledge  of  source  water  character- 
istics (from  in  situ  monitoring  or  historical  data)  ,  stage 
of  development  (inferred  from  satellite  images  and  in  situ 
monitoring),  and  dynamic  processes  (wind  mixing,  advection, 
and  heat  transfer)  a  forecast  of  nutrient  distributions  with 
a  surface  thermal  feature  could  be  made. 
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APPENDIX  A 

Listing   of   Surface   Data:      Time,    Latitude,    Longitude, 
Elapsed   Distance,    Nitrate,    Phosphate,    Nutrient   Ratio, 
ATP,    AATP/ATP,    Chlorophyll,    Temperature 
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APPENDIX   B 

Listing   of  Nansen   Cast   Data:      Station,    Time,    Latitude 
Longitude,    Depth,    Salinity,    Temperature,    Density, 

Nitrate,    ATP 
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APPENDIX   C 
Computer   Programs 
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